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Summary 

The Johnson Space Center performed a prelim- 
inary reduction of the data taken by a Langley- 
designed velocity-measuring system using a 
Hall-effect device during the Crew and Equipment 
Translation Aid (CETA) experiment on the STS-37 
Shuttle flight in April 1991. The data reduction pro- 
vided average velocity profiles. This paper describes 
a new enhanced data reduction technique based on 
a least-squares method that allows instantaneous ve- 
locity profiles to be obtained. Sample velocity pro- 
files of several CETA carts tested in this experiment 
were obtained using this technique. The results show 
that the technique dramatically enhances the veloc- 
ity measurement, especially at the beginning and end 
of the runs. 

Introduction 

A means of moving personnel and equipment on 
the Space Station was investigated on the Crew and 
Equipment Translation Aid (CETA) experiment on 
the STS-37 Shuttle flight in April 1991. Two types of 
instrumentation were used to sense the velocity of a 
cart moving along a specially constructed track. One 
sensor was a Langley-designed special system using 
a Hall-effect device as the basis for a noncontacting 
motion monitor (ref. 1). This system provided the ex- 
act times that the cart passed small magnets placed 
along the track. A reliable average velocity profile 
could also be determined with a minimum amount of 
data manipulation, but with a limited ability to sense 
rapidly changing velocities. Therefore, the velocity 
profiles did not closely follow the true cart velocity, 
especially during cart acceleration and deceleration. 

The other sensor was a three- axis accelerometer 
package installed on the cart. In theory, accelerom- 
eter data can be integrated to provide instantaneous 
velocity profiles at the data system sampling rate. 
The time resolution is limited only by the data sys- 
tem sampling rate and thus can be much higher than 
that obtained from the Hall-effect sensor data (can 
be as high as 300 to 1). However, this method re- 
quires accurate accelerometers and data acquisition 
systems. Information on the initial conditions (ve- 
locity and displacement) and the accelerometer cali- 
bration constants (bias and sensitivity) also need to 
be carefully determined. 

This technique can also be applied in other mea- 
surement areas. For instance, a similar technique can 
be used to verify the performance of the accelerom- 
eters currently being used for critical microgravity 
measurements. Also, the technique can enhance the 
measurements by verifying the calibration constants 
of the instrument in situ. 


Background and Requirement 

The CETA experiment involves accelerating and 
decelerating a cart and its payload along a 50- ft 
track (fig. 1) to various speeds by different propulsion 
methods. The velocity of the cart must be measured 
to evaluate the different propulsion methods as well 
as the effectiveness of the brakes. Because the cart 
does not remain in continuous contact with the track 
in the low-gravity environment in orbit, a noncon- 
tacting type of sensor was required to avoid causing 
drag on the cart, thus contaminating the experiment. 



The sensing system had to be small and battery 
operated (5 V at 50 mA). Another obvious require- 
ment was that the measuring system must be safe 
and not interfere with other subsystems on the Shut- 
tle. Accuracy requirements were that the system 
would provide velocity information up to 6 ft/sec 
within 5 percent of full scale. The output of the sys- 
tem had to be compatible with the battery-operated 
Portable Data Acquisition Package (PDAP) that 
logged all the data on the cart. A means of determin- 
ing the direction of motion was also required (ref. 1). 
The system had to sense the cart motion without 
being affected by the off-axis motion. 

Instrumentation 

The first measurement technique provided used 
a Hall-effect sensor (a solid-state device that senses 
magnetic fields) mounted on the cart to sense mag- 
nets that are mounted in evenly spaced (1-ft) loca- 
tions along the track (fig. 2). Position information 
was recorded directly, and the average velocity was 
determined by measuring the time between encoun- 
ters with consecutive magnets. The Hall-effect sensor 
and the magnets were commercially available. The 
system required no external circuitry and needed no 
adjustments. A second Hall-effect sensor was used 
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Figure 2. Magnet and Hall-probe locations on CETA assembly. 


in tandem to sense the direction of the cart motion. 
This was not essential to the CETA experiment but 
facilitated the data reduction without having to cor- 
relate the acquired data with the video images on 
tape. It also provided a redundant measurement that 
enhanced the system reliability. The second measure- 
ment technique employed three commercially avail- 
able accelerometers to sense the instantaneous accel- 
eration levels of the cart in three axes (longitudinal, 
lateral, and vertical). 

Flight Summary 

The CETA experiment was successfully con- 
ducted on April 8, 1991. Three carts, each using a 
different propulsion method, were successfully oper- 
ated at speeds up to 4 ft /sec. Data were successfully 
acquired on the two Hall sensors and the three ac- 
celerometers. One cart was propelled by manually 
pulling it along the handrail. This method produced 
jerky acceleration but excellent coasting because the 
cart was not in continuous contact with the rail. The 
other carts (mechanical and electrical) had a some- 
what jerky motion and limited coasting because the 
pumping action used a wheel pressed against the 
track. The deceleration in each case was generated 
by a caliper brake controlled by the astronaut. In 
all cases, the maximum acceleration was generated 
during the braking phase of the run. Maximum de- 
celeration rates were approximately 4 ft/sec/sec. 

Data Acquisition 

The raw data of the CETA experiment were 
recorded by the PDAP, which recorded 31 analog 
data channels. The data were sampled 150 times 
per second, digitized to a 12-bit resolution, and time 


tagged. A 40-Hz low-pass Jiher was used to reduce 
the noise and prevent aliasing. A detailed descrip- 
tion of the data acquisition procedure can be found 
in reference 2. - 

Preliminary Data Reduction 

The data of the experiment were initially re- 
duced into two formats: the quick-look and the close- 
up. The gujck-l ook format co ntains data in 1800- sec 
durations over the entire course of the experiment 
(ref. 2) and provides quick overviews of the events 
that happened during the experiment. The time res- 
olution is therefore limited. The close-up format con- 
tains data in 45-sec durations and provides enlarged 
views of the experiment, thus allowing more detailed 
inspections of the data. The close-up data are avail- 
able from the Johnson Space Center (JSC) upon re- 
quest. Figures 3 to 5 (taken from ref. 3), which are 
examples of the close-up data, contain both Hall- 
effect sensor data and accelerometer data obtained 
during the manual cart, mechanical cart, and elec- 
trical cart evaluations. The average velocity profiles 
reduced from these data are also included in these 
figures. 

The average cart velocity is determined by di- 
viding the distance between magnets by the time 
between the Hall-effect sensor transitions. The ac- 
celerometer data can be integrated to produce ve- 
locity profiles. Errors in the~data system and the 
accelerometer are two factors affecting the accuracy 
of the results. The 12-bit data system causes a 
quantization error of about 0.024 percent. Because 
the longitudinal acceleration range of the experi- 
ment is about 0.5# (where 1 g ^ 32.174 ft/sec 2 ), 
the corresponding error is thus 0.00012#. This 
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(b) Reduced forward velocity data. 

Figure 4. Typical Hall-probe and velocity close-up data of channel 29. Data are taken from reference 3. 
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Figure 5. Typical close-up data of accelerometers.. Data are taken from reference 3. 


error can cause a velocity error of 0.08 ft/sec 
(i.e., 0.000120 x 32.2 ft/sec 2 x 20 sec) by the end of 
a typical 20-sec run. The estimated accelerometer 
error is 0.01$, which can cause a 6. 4- ft/sec error by 
the end of a similar run. The maximum root-mean- 
square velocity error due to the data system and ac- 
celerometer is thus 6.4 ft /sec by the end of a 20- sec 
run. 

Enhanced Data Reduction Technique 

This paper presents an enhanced data reduc- 
tion technique that combines the advantages of both 
methods mentioned previously. The data used for 
this technique include the time data, the Hall-effect 
sensor data (channel 29), the velocity data calculated 
from the Hall-effect sensors, and the longitudinal ac- 
celeration data (channel 30). The basic idea is to__ 
use a least-squares method to fit the displacements 
indicated by the Hall-effect sensor data to the cor- 
responding displacements obtained from the double- 
integrated acceleration data. This method provides a 


better estimate of the velocity and displacement pro- 
files. The resulting estimated parameters include thj 
initial conditions of the velocity and displacement, 
and also the accelerometer sensitivity and bias. A 
detailed derivation of the equations involved in this 
method is provided in the appendix. A description 
of the procedure for this enhanced data reduction 
technique is given in the following six steps: ; 

1. Locate the time window from the close-up 
data. 

The time window desired is first located in graphs 
designated by “channel 28 aft speed” and “channel 29 
forward speed” of the close-up data (similar to those 
in figs. 3 and 4). These graphs contain the velocity 
profiles obtained from reducing the two sets of Hall- 
effect s enso r data. Because the starting point of 
each test is not obvious from the data, the desired 
starting point of the time window is assumed to be 
approximately 1 sec before the beginning of each 
velocity profile. Data in the opposite direction from 
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Figure 6. Forward velocity profile of manual cart. Constants at to — 3120 sec: so 0.56 ft, t>o 0.87 ft/sec; b 0.27 ft/sec ; 
m = -30.82 ft/sec 2 /# units. 



Figure 7. Forward velocity profile of mechanical cart. Constants at to = 4720 sec: so — 0.33 ft; uo — 0.25 ft/sec; b 0.30 ft/sec , 

m = -31.18 ft/sec 2 /g units. 



Figure 8. Forward velocity profile of electrical cart. Constants at to = 7004 sec: sq — -0.73 ft; vq - 0.29 ft/sec; b 0.28 ft/sec ; 

m = -29.00 ft/sec 2 /p units. 
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the previous run should never be included because 
these data produce anomalous results. 

2. Extract the time data, Hall-effect sensor data, 
and acceleration data. 

The time data, the Hall-effect sensor data 
recorded from channel 28, and the longitudinal ac- 
celeration data of the corresponding time window 
are extracted from the standard data file (SDF) and 
stored as a single data file. The SDF is the principal 
processed file that contains the CETA experiment 
data. (See ref. 3.) 

3. Import the data file into spreadsheet software. 

4. Execute the macro within the spreadsheet. 

The spreadsheet contains a macro (a set of com- 
mands that can be executed with a single command) 
that integrates and then extracts the displacements 
from the acceleration data. The macro first numer- 
ically integrates the acceleration data to obtain the 
velocity and displacement profiles (zero initial condi- 
tions). The gradients of the Hall-effect sensor data 
are also calculated. The gradients of the data^re 
then sorted for magnitudes over 100 V/sec. This 
magnitude is chosen as a criterion to detect the tran- 
sition periods in the data. (Transition periods are 
the time intervals during which the Hall-effect sen- 
sor changes state as it passes over a magnet.) The 
macro then looks for sign changes in the computed 
gradients. Whenever a sign change occurs, the macro 
extracts the corresponding time and displacements 
from the integrated acceleration data. 

5. Calculate the matrix elements and the curve- 
fitting constants. 

The data extracted in step 4 are then used to 
least-squares fit the displacements indicated by the 
Hall-effect sensor, which are relative to the first tran- 
sition. The distance between two consecutive tran- 
sitions is 1 ft. The least-squares fit is done by first 
calculating the matrix elements in equation (A16) 
shown in the appendix. Equation (Al7) is then eval- 
uated to obtain the initial parameters sq, uq, 6, and 
m, where so is the initial displacement relative to the 
first magnet encountered by the vehicle, vq is thejni- 
tial velocity of the time window, and b and m are the 
accelerometer bias and sensitivity correction factors, 
respectively. 

6. Reintegrate the acceleration data to obtain the 
instantaneous velocity and displacement profiles. 

The acceleration data are reintegrated using equa- 
tions (A6) and (A7) in the appendix with the param- 


eters obtained in step 5 to correct the velocity and 
displacement profiles. 


Results 

Figures 6, 7, and 8 illustrate the velocity pro- 
files obtained by using the enhanced data reduction 
technique for the manual, mechanical, and electri- 
cal carts, respectively, tested in the CETA experi- 
ment. Shown in each figure is the corrected velocity 
profile using the enhanced data reduction technique 
compared with the velocity profile obtained from the 
Hall-effect sensor data. The parameters (sq, t>o, 6, 
and m) estimated by this least-squares method are 
also listed. Because the acceleration data are used, 
the instantaneous velocities are now available instead 
of the average velocities. Also, th e enhanced resolu- 
tion due to using the acceleration data provides more 
detailed velocity information in the beginning, end, 
and between magnets. The error contributed by this 
technique is primarily from the numerical integra- 
tions of the discrete acceleration data. 


Concluding Remarks 

The newly developed enhanced data reduction 
technique provides an improved procedure that al- 
lows least-squares minimization to become possible 
between data sets with an unequal number of data 
points. This technique was applied in the Crew and 
Equipment Translation Aid (CETA) experiment on 
the STS-37 Shuttle flight in April 1991 to obtain the 
velocity profile from the acceleration data. This new 
technique uses a least-squares method to estimate the 
initial conditions and calibration constants. These 
initial conditions are estimated by least-squares fit- 
ting the displacements indicated by the Hall-effect 
sensor data to the corresponding displacements ob- 
tained from integrating the acceleration data. The 
velocity and displacement profiles can then be re- 
calculated from the corresponding acceleration data 
using the estimated parameters. The technique en- 
ables instantaneous velocities to be obtained from the 
test data instead of only average velocities at vary- 
ing discrete times. Therefore, more detailed velocity 
information is offered, particularly during periods of 
large acceleration or deceleration. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
March 31, 1993 
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Appendix 

Derivation of the Least-Squares-Fit 
Matrix Equation 

The new enhanced data reduction technique is 
designed to obtain least-squares estimates of the cor- 
rection factors and initial conditions of the accelera- 
tion calibration by fitting the integrated acceleration 
data to the Hall-effect sensor displacement data. 

Let x(t) denote the uncorrected acceleration data 
obtained as a function of time during a CETA run. If 
the calibration factors (the sensitivity m and the bias 
b ) are in error, the corrected acceleration data (a(t)) 
may be obtained from x(t) by correction factors as 
shown in the equation 

a(t ) = m x(t) + b (Al) 

The corrected velocity data {v(t)) are obtained by 
integrating the corrected acceleration data to give 

v(t) = uo + f a{t) dt (A2) 

Jto 

where vq and to denote the initial velocity and time, 
respectively. 

Substituting equation (Al) into equation (A2) 
gives 

v(t) — vo + [ [m x(t) 4- b]dt (A3) 

Jt 0 

The result of the integration is the equation 

v(t ) = vq + b(t -to) + m /'(t) (A4) 

where 



For discrete data, equations (A4) and (A5) can 
be evaluated at sampling time ( t n ) as, respectively, 

v(tfi) = vo + b{tn — to) th I (tn) (A6) 

and 

s(t n ) = So-f^o^n+to)"^- -b(t n + to) 2 + m I(t n ) (A7) 
where 

t n = nT (n = 0, 1 , N ) 

T denotes the sampling period of the data system, 
and N denotes the number of data points. The 
terms I l (t n ) and J(t n ) are obtained by numerical 
integration. In this paper, rectangular integration is 
used. The inner integral (J'(t)) is thus approximated 

by 

k 

l'(t k ) « ^2 At* x(U) {k = 1, 2, N) (A8) 

i = 1 

where A t^ = tj — t^_i and /'(to) = 0. The term J(t n ) 
is approximated by rectangular integration as 

k k j 

I(t k ) « £ A tj I'(tj ) A h At i *(<<) ( A9 ) 

j= 1 J =1 *=1 

Let s h {t H j) denote the displacement indicated by the 
Hall-effect sensor when it passes a magnet at time 
t Hy The set of Hall-effect sensor sampling times 
(i.e., where j = 1, 2, ..., M) is a subset of the 

data system sampling times (i.e., {t*} where i = 1, 
2, N). The sum of squares of differences between 

the corrected displacement data and the data from 
the Hall-effect sensor, summed over the Hall-effect 
sensor sampling times, is described by the equation 


Equation (A4) can be further integrated to obtain 
the corrected displacement data ( s(t )). That is, 


s(t) = sq + f v(t) dt 

Jt 0 


= s 0 + / [vo + b(t - to)]dt + m 

Jt 0 

= sq + uo(* - to) + “ *o) 2 + m ^(0 (AS) 



M 

e 2 = £>(*//,•) - s h {t H j)? (M < N ) (A10) 

j= 1 

where M denotes the number of Hall-effect sensor 
samples. 

If we substitute equation (A7) into equa- 
tion (A10) and let tjjj - to be replaced by At^, 
the equation then becomes 

M l 

£? = ^ ^ [so 4- t>o At//j 4 — b At 4 tti Z(t//j) — s h(t//j)] 
j=\ 

(AH) 


where so denotes the initial displacement and 

I(t) = f I'(t) dt 
Jto 
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The least-squares estimates of parameters sq, uq, b, and m are obtained by minimizing e 2 . That is, 


a £= a 

(A12) 


(A13) 

o 

II 

cs 1.0 

(A14) 

¥- =o 

(A15) 


If we substitute equation (All) into equations (A12), (A13), (A14), and (A15), the corresponding partial 
derivatives are, respectively, 

M j 

2^ [so + vo At xj + -b At jjj + m ~ — 0 

i-} 

M L 

2 5Z ( s 0 + «0 At Ifj -I- -b At 2 Hj + m T(t H j) - Sfc(t^j)] = 0 
3 - 1 

M l 

2 yz l s 0 + VQ At H j + -b Atffj + m I(t H j) - s h (t H j)]At H j = 0 
i= 1 


M 


2 J2 ( s 0 + vq At Hj + -b At 2 Hj + m I(t Hj ) - s h (t Hj )] 2 I(t H j) = 0 

j=l 


The matrix form of this set of equations is given by 


N 

M 

E At Hj 

Af 

E 

J=l 

M 1 

E 

j=l 



r m i 

E s ^(t//j) 

i=i 

Af 

E At Hj 

j=l 

M 

E At lj 

E At », 

J=1 

Af 

E Atjyj- I(ttfj) 
i=i 


“so" 

VO 


E At Hj s /i(tf/j) 
i=i 

M 

E At h 

j = i 

Af 

E 

J = 1 

E Atjy ■ 

J=1 

Af 

E At/fj I(tffj) 
j=l 


b 

_m. 


Af 

E At ff ^ s/j(t//j) 

j=l 

M 

E 

Li=! 

M 

E Attfj I(t^) 

J=1 

M n 

E At^ 

i=i 

1 

-4-» 

CN 

l-H 



Af 

E sft(t/fj) 

. j = i J 


(A16) 


8 



The parameters can thus be obtained by manipulating the matrix into the form 


N 

At 

E At Hj 

3= 1 

At 

At 

E 

j=i 

-1 

r M 1 

E sfrOtfj) 
i=i 

At 

E At Hj 

i=i 

<1 

E 

3= 1 

At 

E At//, I(t//j) 
i=i 


M 

J=1 

At 
3 = 1 

M 

E At^. 

i=i 

M 

E At 4 H ■ 
i=i 

At 

E At H j I(t H j) 

j= i 


Af 

E s /i ("t jr ) 

j=l 

M 

E i(t^) 

At 

E Atf/j I(t//j) 

J=1 

A/ 

E At^ I(tffj) 

i=i 

At 

E 

i=i 


U 

L j=i -1 

/ A 


(A17) 


The parameters thus estimated can then be applied to equations (A6) and (A7) to obtain corrected velocity 
and displacement profiles. 
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